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Abstract 

The performance of multiple-input multiple-output wireless systems is investigated in the presence of 
statistical queueing constraints. Queuing constraints are imposed as limitations on buffer violation probabilities. 
The performance under such constraints is captured through the effective capacity formulation. A detailed analysis 
of the effective capacity is carried out in the low-power, wideband, and high-signal-to-noise ratio (SNR) regimes. 
In the low-power analysis, expressions for the first and second derivatives of the effective capacity with respect 
to SNR at SNR = are obtained under various assumptions on the degree of channel state information at the 
transmitter. Transmission strategies that are optimal in the sense of achieving the first and second derivatives 
are identified. It is shown that while the first derivative does not get affected by the presence of queueing 
constraints, the second derivative gets smaller as the constraints become more stringent. Through the energy 
efficiency analysis, this is shown to imply that the minimum bit energy requirements do not change with more 
strict limitations but the wideband slope diminishes. Similar results are obtained in the wideband regime if 
rich multipath fading is being experienced. On the other hand, sparse multipath fading with bounded number 
of degrees of freedom is shown to increase the minimum bit energy requirements in the presence of queueing 
constraints. Following the low-SNR study, the impact of buffer limitations on the high-SNR performance is 
quantified by analyzing the high-SNR slope and the power offset in Rayleigh fading channels. Finally, numerical 
results are provided to illustrate the theoretical findings, and to demonstrate the interactions between the queueing 
constraints and spatial dimensions over a wide range of SNR values. 

I. Introduction 

Having multiple antennas at the transmitter and receiver has been shown to improve the performance 
significantly in terms of both reliability and throughput when the channel fading coefficients are known 
at the receiver and/or transmitter. Due to these promising gains in the performance, information-theoretic 
analysis of multiple-input multiple-output (MIMO) channels has attracted much interest in the research 
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community. In particular, considerable effort has been expended in the study of the capacity of MIMO 
wireless channels (see e.g., [1] and the references therein). For instance, multiple-antenna capacity is 
studied in the low-power regime in [2] and [3], and in the high-SNR regime in [4]. In most studies on 
MIMO channel capacity, ergodic Shannon capacity formulation is employed as the main performance 
metric. However, this formulation does not capture the performance in the presence of quality-of-service 
(QoS) limitations in the form of constraints on queueing delays or queue lengths, although providing 
QoS assurances is of paramount importance in many delay-sensitive wireless systems, e.g., voice over 
IP (VoIP), and interactive and streaming video applications. 

In [5], effective capacity is proposed as a metric that can be employed to measure the performance 
in the presence of statistical QoS limitations. Effective capacity formulation uses the large deviations 
theory and incorporates the statistical QoS constraints by capturing the rate of decay of the buffer 
occupancy probability for large queue lengths. Hence, effective capacity can be regarded as the maximum 
throughput of a system operating under limitations on the buffer violation probability. This formulation 
is tightly linked and in a sense dual to the concept of effective bandwidth [6] [7] that is employed in 
the analysis of how much resource in terms of service rates is needed to support a given time-varying 
arrival process. The analysis of the effective capacity in various wireless communication settings has 
been conducted in several recent studies (see e.g., [9] - [16]). 

In this paper, we study the effective capacity of MIMO wireless channels. In particular, we consider 
the low-power, wideband, and high-SNR regimes and identify the impact of the QoS limitations^ on 
the performance. We would like to note that recently references [17] and [18] have also investigated 
the effective capacity of multiple- antenna channels. In [17], the authors study the multiple-input single- 
output (MISO) channels and determine the optimal transmit strategies with co variance feedback. In 
[18], the concentration is on the MISO and single-input multiple-output (SEVIO) channels. Analysis of 
MIMO channels is carried out only in the large antenna regime in which the number of receive and/or 
transmit antennas increase without bound. In addition, the authors in [18] consider a MIMO channel 
matrix with independent and identically distributed (i.i.d.) zero-mean Gaussian entries, and consider 
equal power allocation across the antennas. In this paper, we consider a general MIMO link model in 
which the fading coefficients have arbitrary distributions and are possibly correlated^, provide a detailed 
study of the low-power, wideband, and high-SNR regimes, investigate the transmission strategies under 
various assumptions on the degree of channel knowledge at the transmitter, and identify the impact of 



'Throughout the paper, we use the terms "QoS constraints", "queueing constraints", and "buffer constraints" interchangeably. 
2 Only in the high-SNR regime, we coi 
unit-variance, Gaussian random variables. 



2 Only in the high-SNR regime, we concentrate on the canonical MIMO model in which the fading coefficients are i.i.d. zero-mean, 
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QoS constraints on the performance. The original contributions of this paper are the following: 

1) We obtain expressions for the first and second derivatives of the effective capacity at snr = 
under various assumptions on the availability of channel knowledge at the transmitter, and show 
that while the first derivative is independent of the queueing constraints, the second derivative 
diminishes as the constraints become more stringent. Transmission strategies that achieve these 
derivatives are identified. 

2) As a result of the findings on the derivatives of the effective capacity, we determine in the 
low-power regime that the minimum bit energy requirements in the presence of QoS limitations 
are the same as those attained in the absence of such constraints. On the other hand, we show 
that the wideband slope decreases under more strict queueing constraints, indicating that energy 
expenditure increases unless one is operating at the minimum bit energy level. 

3) Under certain assumptions, we show that the results obtained in the low-power regime apply to 
the wideband regime with rich multipath fading. In contrast, we establish that sparse multipath 
fading has a significant impact on the performance in the wideband regime. In particular, we 
prove that minimum bit energies greater than that achieved in the absence of QoS constraints are 
required if the number of degrees of freedom in the form of noninteracting subchannels remain 
bounded as the bandwidth increases. 

4) Considering i.i.d. Rayleigh fading channel model, we identify the effect of QoS limitations on the 
performance in the high-SNR regime by determining the high-SNR slope and power offset values. 

The organization of the rest of the paper is as follows. We describe the MIMO channel model in 
Section UH In Section HHl we provide a description of the effective capacity formulation, and apply 
it to the MIMO setting. In Section [IV] we study the effective capacity in the low-power regime and 
determine the first and second derivatives of the effective capacity at zero snr. Subsequently, we apply 
the derivative expressions to investigate the energy efficiency. In Section |V1 we explore the effect of QoS 
limitations in the wideband regime, and identify the minimum bit energy requirements. In Section [VH 
we concentrate on the high-SNR regime, and determine the impact of QoS constraints on the performance 
in the i.i.d. Rayleigh fading channel. Finally, we provide numerical results in Section [VTll and conclude 
in Section [VTTT1 

II. Channel Model 

We consider a MIMO channel model and assume that the transmitter and receiver are equipped with 
n T and n R antennas, respectively. Assuming flat-fading, we can express the channel input-output relation 
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as 

y = Hx + n. (1) 

Above, x denotes the ut x 1-dimensional transmitted signal vector, and y represents the ur x 1- 
dimensional received signal vector. The channel input is assumed to be subject to the following average 
energy constraint: 

E{||x|| 2 }<§ (2) 

where B is the bandwidth of the system. When the bandwidth is B, we can assume that B input vectors 
are transmitted every second, and © implies that the average power of the system is limited by P. In 
(QQ), n with dimension n R x 1 is a zero-mean Gaussian random vector with E{nn^} = N I, where I is 
the identity matrix. The signal-to-noise ratio is defined as 

_ E{||x|| 2 } P 
SNR "E{||nP} "r^AV (3) 

We also define the normalized input covariance matrix as 

K = (4) 
P/B ' K) 

Note that the average energy constraint in © implies that the trace of the normalized covariance matrix 
is upper bounded by 

tv{K x )<l. (5) 

Finally, in CD, H denotes the n R x n T -dimensional random channel matrix whose components are 
the fading coefficients between the corresponding antennas at the transmitting and receiving ends. 
Unless specified otherwise, the components of H are assumed to have arbitrary distributions with finite 
variances. Additionally, we consider the block-fading scenario and assume that the realization of the 
matrix H remains fixed over a block of duration T seconds and changes independently from one block 
to another. 

III. Effective Capacity of a MIMO Link 

In [5], Wu and Negi defined the effective capacity as the maximum constant arrival rate that a given 
service process can support in order to guarantee a statistical QoS requirement specified by the QoS 



4 



exponent 9 If we define Q as the stationary queue length, then 9 is the decay rate of the tail of the 
distribution of the queue length Q: 

Um l°gP(Q>9) = _ 

q->oo q 

Therefore, for large q max , we have the following approximation for the buffer violation probability: 
P(Q > Qmax) ~ e~ e?max . Hence, while larger 9 corresponds to more strict QoS constraints, smaller 
9 implies looser QoS guarantees. Similarly, if D denotes the steady-state delay experienced in the 
buffer, then P(D > d max ) « e~ eSdmax for large d max , where 5 is determined by the arrival and service 
processes [11]. Therefore, effective capacity formulation provides the maximum constant arrival rates 
that can be supported by the time-varying wireless channel under the queue length constraint P(Q > 
Qmax) < e~ eqmax for large q max or the delay constraint P(D > d max ) < e ~ 9Sdraa * for large d max . Since 
the average arrival rate is equal to the average departure rate when the queue is in steady-state [8], 
effective capacity can also be seen as the maximum throughput in the presence of such constraints. 
The effective capacity is given by ([5], [6], [7]) 

A{ ~ 6) - limllog e E{e- e5 W} (7) 



9 t^oc 9t 

where S[t] = ^2 l i=1 R[i] is the time-accumulated service process and {_R[z],z = 1,2,...} denotes the 
discrete-time stationary and ergodic stochastic service process. Under the block-fading assumption, the 
effective capacity formulation simplifies to 

AM) - ^log e E{e-^}. (8) 



9 9T 

Under a short-term power constraint, the stochastic service process in a MIMO channel with a given 
normalized input covariance matrix is 



B log 2 det ( I + -^-HK 3 .H t ] = B log 2 det (i + ra^SNRHK x H^) bits/s 



(9) 



where B denotes the bandwidth of the system and snr is as defined in ©. We first consider the case in 
which H is perfectly-known at the receiver and transmitter. In this scenario, the transmitter can adapt 
the input covariance matrix with respect to each realization of H in order to maximize the service 
rate. Therefore, using the formulation in ([8]), we can express the effective capacity normalized by the 

3 For time-varying arrival rates, effective capacity specifies the effective bandwidth of the arrival process that can be supported by the 
channel. 
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bandwidth and the receive dimensions as 

C s (snr, 9) = — — 4 — log.E { exp I — 9TB max log 2 det (I + to^snrHK^H 1 ") I \ bits/s/Hz/dimension 

a(K x )<i 

(10) 

As 9 vanishes, the QoS constraints become loose and it can be easily verified that the effective capacity 
approaches the ergodic channel capacity, i.e., 

lim C £ (snr, 9) = — E I max log 2 det (I + to^snrHK^H*) } . (11) 

'tr(K~)<l 

For 9 > 0, the effective capacity is in general smaller than the ergodic capacity. We can easily see this 
by interchanging the logarithm and the expectation in (flOl ) and applying the Jensen's inequality: 

C b (snr, 9) = - 6T l BnR log e E <j exp | -9TB max log 2 det (i + torSNrHK,!^) ] ) (12) 
< - QT l Bn E < log e exp I -9TB max log 2 det (i + t^snrHK^H 1 ") ] ) (13) 

R I \ tr(KaO<l 

= — E < max log 2 det (I + to^snrHK^H 1 ") \ . (14) 

TOr K x ^0 1 

^tr(K x )<l 

Above, we have assumed that H is perfectly known at the transmitter. If, on the other hand, only 
statistical information regarding H is available at the transmitter, then the input covariance matrix can 
be chosen to maximize the effective capacity. In such a case, the normalized effective capacity can be 
expressed as 

C b (snr, 9) = max -77=-^ — l°g e E (exp (-9TB log 2 det (i + to^snrHK^H 1 ")) } bits/s/Hz/dimension. 

k^o 9TBriji 

tr(K x )<l 

(15) 

For a given (and not necessarily optimal) input covariance matrix K. x , we call the throughput as effective 
rate and express it as 

Re(snr, 9) = log e E {exp {-9TB log 2 det (I + to^snrHK^H 1 ")) } bits/s/Hz/dimension. 

9TBtir 1 v v ' ' ' 

(16) 

In practice, uniform power allocation across the antennas might be preferred. In this case, K. x = ^1, 
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and the effective rate can be written as 

R E idfsNR, 6) = -—— log e E | exp ( -9TB log 2 det (I + — ^snrHHM ) 1 bits/s/Hz/dimension 

0TBn R y V " V n T J J ) 

(17) 

where the subscript "id" is introduced to denote that this expression is the throughput when the 
covariance matrix is proportional to an identity matrix. 

Note that the effective capacity and effective rate expressions in (flOl) . (fT5l) . (fT6l) . and (fTTT) are 
proportional to the logarithm of the moment generating function of the instantaneous transmission 
rates. 

Since the subsequent analysis assumes that the QoS exponent is fixed as power diminishes or increases 
or bandwidth increases, we generally suppress the argument 9 and write the effective capacity and rate 
as C^(snr) and Re(snr), respectively. 

Finally, before we go through a more detailed analysis of the effective capacity in the following 
sections, we would like to discuss several implicit assumptions made in the formulations provided in 
this section. The service rate expression in © implies that the maximum transmission rates are equal 
to the instantaneous channel capacity in each block of duration T. Hence, we implicitly assume that the 
number of symbols in each block, TB, is large enough for this assumption to have operational meaning 
in practice. In (fT5l) . it is assumed that the service rate is still given by © and hence the transmitter 
employs variable-rate transmission scheme, even though the transmitter does not know the instantaneous 
realizations of H. Note this can be accomplished by using recently developed rateless codes such as 
LT [19] or Raptor [20] codes, which enable the transmitter to adapt its rate to the channel realization 
without requiring CSI at the transmitter side [21], [22]. It is also important to note that the analysis 
conducted in this paper apply in the large-queue-length regime. If the buffer size is finite and small, 
then the arrival rates that can be supported by the system will be smaller than those considered in the 
paper, and in this case, one has to consider packet loss probabilities as well. Therefore, if the above- 
mentioned conditions and assumptions are not satisfied in the system, then the performance degradation 
will be more severe. For such cases, the results of this paper can be seen as fundamental limits (or 
upper bounds) which can serve as benchmarks for system performance. 
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IV. Effective Capacity in the Low-Power Regime 

A. First and Second Derivatives of the Effective Capacity 

In this section, we study the effective capacity in the Iow-snr regime and investigate the impact of 
the QoS exponent 8. In particular, we consider the following second-order expansion of the effective 
capacity under different assumptions on the degree of channel state information: 

SNR^ 

C e (snr) = C b (0)snr + C E (0)— - + o(snr 2 ) (18) 

where Qe(0) and Qe(0) denote the first and second derivatives of the effective capacity with respect to 
snr at snr = 0. We first have the following result when the channel is perfectly known at the transmitter 
and receiver. 

Theorem, 1: Assume that the realizations of the channel matrix H are perfectly known at the receiver 
and transmitter. Assume further that the transmitter is subject to a short-term power constraint and 
hence is not allowed to perform power adaptation over time. Then, the first and second derivatives of 
the effective capacity in (flOl) with respect to snr at snr = are 

Ce(0) = -^-E{A max (HtH)} (19) 

log e 2 

and 

C*(0) = 6 -^f [E 2 {A max (HtH)} - E{A 2 nax (H t H)}] - tt— ^-xE{A 2 nax (H t H)} (20) 
log e 2 flog e 2 

where A max (H^H) denotes the maximum eigenvalue of H^H, and I is the multiplicity of A max (H^H). 
Proof: For a given input covariance matrix K. x , the effective rate is expressed as 

R b (snr) = -—— log e E {exp {-8TB log 2 det (I + n^SNRHK^Ht) ) } (21) 

81 BriR 

1 log e E{exp(-0T£log 2 det(I + n K SNR$))} (22) 



8TBn R 



8TmT R loge E { 6XP [~ 9TB S log 2 + nRSNRX^))] | (23) 



1 f I 8TB 

log e E < exp I !og e (1 + n fl SNRA 4 ($)) ) \ (24) 



8TBn R ° e 1 "\ log e 2 
1 



8TBn R 



log e E{/(sNR,0)}. (25) 
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In ([221 above, we have defined $ = HK x Ht ([231) is obtained by noting that the determinant of a 
matrix is equal to the product of its eigenvalues, i.e., det (I + n#SNR<I>) = Yl^l + n R SNR\i(Q>)), and 
also using the fact that the logarithm of a product is equal to the sum of the logarithms of the terms 
in the product. In (|24)) . the base of the logarithm is changed from 2 to e. In (|25l) . we have defined the 
function /(snr, 9) = exp f-g^ !og e (1 + n;?SNRAi($))J. 
Now, taking the derivative of R E with respect to snr yields 

w - - ^n R E{/(sNR, 9) } E {"S^ i + ££i ( *) /(sNM) } • <26) 

Noting that the function / evaluated at snr = is one, i.e., f(0,9) = 1, we can easily see from (|26|) 
that the value of the first derivative of the effective rate at snr = is 

Mo) = ^e{x>(*)} = j^M*)} = To^ E { tr ( HK ^ Ht )} W 

where we have used the fact that the sum of the eigenvalues of a matrix is equal to its trace. Note that 
the normalized input covariance matrix is by definition a positive semidefinite Hermitian matrix. 
As a Hermitian matrix, K x can be written as [31, Theorem 4.1.5] 

K^UAU^^g^uJ (28) 

i=i 

where U is a unitary matrix, {uj} are the column vectors of U and form an orthonormal set, A is 
a real diagonal matrix, {di} are the diagonal components of A. Since K x is positive semidefinite, we 
have di > 0. Moreover, since all available energy should be used for transmission (i.e., the average 
energy and hence trace constraints should be satisfied with equality), we have tr (K^) = Yll=i d% = 1- 
Combining (|27b and (128b . we can now write 

~ ' ' (29) 



km = io^2 E ^ r(HK;EHt) ^ = io^E^ E { tr ( Hu ^ tHt )} 



i=i 



log, 



— ^diEjujHtHin} (30) 



e ~ i=i 



< E{A max (H t H)}. (31) 
log e 2 

where A max (H^H) denotes the maximum eigenvalue of the matrix H^H. The upper bound in (l3Tb 
follows from the facts that di £ [0, 1] and Ei<^ = 1> an d from [31, Theorem 4.2.2] which states that 
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since H+H is a Hermitian matrix and {u;} are unit vectors, we have 

ujH+Hu* < A max (H+H) Vi. (32) 

The upper bound in (I3TI) can be achieved by beamforming in the direction in which A max (H^H) is 
achieved, i.e., by choosing the normalized input covariance matrix as 

= uu f (33) 

where u is the unit-norm eigenvector that corresponds to the maximum eigenvalue A max (H^H). This 
lets us conclude that 

Cs(0) = r ^-E{A max (HtH)} (34) 
log e 2 

proving (fT9l) . 

Before proceeding to the proof of the second derivative result, we would like to note that transmission 
in the maximal-eigenvalue eigenspace of H^H is indeed necessary to achieve the first derivative. 
Therefore, it is also necessary to attain the second derivative of the effective capacity at zero snr. 
In a general scenario in which A max (H^H) has a multiplicity of I > 1, an input covariance matrix in 
the following form is required: 

l 

K x = ^2 a i u i u l (35) 

i=i 

where aj E [0, 1] and Y!h=i a i = 1' an( ^ { u «}i=i are the orthonormal eigenvectors that span the maximal- 
eigenvalue eigenspace of H^H. 

Now, we turn to the analysis of the second derivative. Differentiating R E in (|26|) once more with 
respect to snr, we obtain 



Re(SNR) = " 5^2 g{/(s»M)} ~ E | £ l +nr SNRA i( $) /(SNR - 9) 

1 1 _ -n,A?(* 



+ i^ E{/( S NM)} E {^ (i + "^A, ( ; )P /(^- ( ')^ ™ 

9TBn R 1 



log* 2 E{/(SNR,0)} 




/(snr, (9) 
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Again noting that f(0,9) = 1, we have 

^ (0) = i$r ( E2 {l>(*>} - E { (sx* 

= ^Y(E 2 {tr(*)}-E{tr 2 ($)})- i ^E{tr($t$)}. 

In obtaining (1381) . we have used the facts that = tr(<&) and = tr($ t $). 

As described above, an input covariance matrix that is in the form given in ([35]) is required to achieve 
the second derivative of the effective capacity at snr = 0. For such a covariance matrix, it can be easily 
verified that 

E{tr($)} = E{tr(HK x .H f )} = E |A max (H t H)} (39) 

and 

E{tr($ f $)} = E {tr (HKj.tfHK^tf) } = E j c^luJl^Huil 2 j (40) 

I i,j ) 

= E |a^(H+H) ^ a^lutu^j (41) 

= E|AL x (H t H)ga l 2 | (42) 

>yE{AL x (HtH)} (43) 

where (|4T|) follows from the fact that {u^} are the eigenvectors that correspond to A max (H^H) and 
hence H*Hu, = A max (H^H)uj, (1421) follows from the orthonormality of {uj} which implies that 

+ 1 if % — j 

u]in ={ . (44) 

[ if i± j 

Finally, (|43l follows from the properties that G [0, 1] and Yl\=i OL % — ^-> an d the fact that Y^i=i a i 
under these properties is minimized by choosing ctj = j, which leads to the lower bound Yl\=\ a i — f • 
We note from (|39l) that given the required the covariance structure in (1351) . the first term in the expres- 
sionof R £ (0) in® is fgf (E 2 {*(<&)} - E{tr 2 ($)}) = fgf (E 2 (A max (HtH)} - E (A 2 nax (HtH)}) 
for all possible On the other hand, the second term in (1381) is minimized by having = j for all i, 
i.e., by equally allocating the power in the orthogonal directions in the maximal-eigenvalue eigenspace. 
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Therefore, the input covariance matrix K x = j Yl\=\ u « u ! maximizes Re(0), and we have 

C B (0) = (E 2 (A max (H t H)} - E (A 2 nax (H t H)}) - {A^H)} (45) 

proving (|20l ). ■ 
Next, we consider the case in which the transmitter has only statistical knowledge of the channel. 
Theorem 2: Assume that while the receiver perfectly knows the channel matrix H, the transmitter 

only has the knowledge of E{H^H}. Then, the first and second derivatives of the effective capacity in 

(fT5l) are 

Cs(0) = -^-A max (E{HtH}) (46) 
log e 2 



and 



C*(0) = ^^ALMHtH}) - mm j^a^ ( 9 ^^E{(u\H^Hu l )(u]^nu J )} + -^-E{|u]Ht H u, 
l0ge2 Vlog e 2 log e 2 

(47) 

where A max (E{H^H}) denotes the maximum eigenvalue of E{H^H}, and / is the multiplicity of 
A max (£{HtH}). 

Proof: Note from (l30l) that for a given covariance matrix K x = UAU^ = XT=i ^i u i u l> me fi rs t 
derivative of the effective rate is 

MO) = d ^ {ulHtHiii} (48) 

i=l 

= ; r V ^ujE {HtH} u, (49) 

lo g e 2 ^ 

< r ^-A max (E{HtH}) (50) 

where (|49l follows by noting that the transmitter has only statistical knowledge of H, and the input 
covariance matrix and hence {iij} cannot depend on the realizations of H. Therefore, {u,} are deter- 
ministic and can be taken out of the expectation. Now, the upper bound in (l50l) . similarly as discussed in 
the proof of Theorem [H is achieved by transmitting in the maximal-eigenvalue eigenspace of E{H^H}. 
Therefore, a covariance matrix in the following form is required to achieve the first derivative of the 
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effective capacity: 

i 

Kj. = OiUiiiJ (51) 

i=i 

where {uj} are the orthonormal eigenvectors spanning the maximal-eigenvalue eigenspace of E{H^H}, 
I is the multiplicity of A max (E{H^H}), and {a,} are constants taking values in [0, 1] and has unit sum, 
i.e., Yh=i a i = 1- Consequently, this covariance structure is also necessary to attain the second derivative 
of the effective capacity. Employing the second derivative expression in (1381) with the covariance matrix 
in (IBTI) . and maximizing Re(0) with respect to all possible choices of we easily obtain (|47T) . ■ 
Using the results seen in the proofs of Theorems Q] and [2l we can also immediately obtain the 
following result when the power is uniformly distributed across the transmit antennas and hence we 
have K. x = —I. 

TIT 

Corollary 1: Assume that the input covariance matrix is K x = ^1. Then, the first and second 
derivatives of the effective rate R^,id given in (fTTT) are 

Rem(P) = — ^— E{tr(HtH)} (52) 
n T log e 2 

and 

Rb#(P) = ^TT% [E 2 {tr(Ht H )} -E{tr 2 (Ht H )}] - — E{tr ((H t H) 2 )}. (53) 
nip logg 2 njAog e 2 

Remark 1: Note that the common theme in the results of Theorems Q] and [2l and Corollary Q] is that the 
first derivative does not depend on 9 and hence does not get affected by the presence of QoS constraints. 
Indeed, the first derivative expressions are equal to the ones obtained when Shannon capacity, rather 
than effective capacity, is considered [2]. On the other hand, the second derivative is a function of 9 and 
in general decreases as 9 increases or equivalently as the queueing constraints become more stringenj^. 

B. Energy Efficiency in the Low-Power Regime 

The expressions of the first and second derivatives enable us to analyze the energy efficiency in the 
low-power regime. The minimum bit energy under QoS constraints is given by [2] 

El = , im _^_ = J_. (54) 

N 0min snr^oC s (snr) C e (0) 

4 Note that E 2 {A max (H t H)} < E{ALx(H t H)} and E 2 {tr(H t H)} < E{tr 2 (H f H)}. 
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At j*- . , the slope S of the spectral efficiency versus E b /N (in dB) curve is defined as [2] 



Nomm 



^ , 1,111 — e^ E( ":1 e, 1»1"H,„2. <55) 



^li^l 10 login # - 10 log ln ^ 



Considering the expression for normalized effective capacity, the wideband slope can be found from 

2(C (0)) 2 

S = 1 . EK >> log e 2 bits/s/Hz/(3 dB)/receive antenna. (56) 

— Ce(0) 



Corollary 2: Applying the results of Theorem Q] to the above formulation, we obtain 

E b log. 2 



iV 0min E{A max (HtH)} 



(57) 



c = 2E 2 {A max (H t H)} 

° T E { A Lx(HtH)} + ff^(E{AL x (HtH)}-EHA max (HtH)}) 



2 



^f/<<w(H)) + ^ (*(<w(H)) 



(59) 



log e 2 

where ft(cr max (if)) is the kurtosis of maximum singular value of the matrix H and is defined as 



_ EKx(H)} E{A 2 nax (HtH)} 
1 maxl jj E 2 {aL x (H)} E 2 {A max (HtH)}' (bU) 
Remark 2: In [2], Shannon capacity is considered and it is shown that ^ = „ f . "^V,*^ and 

1 J JVo m i n JtL{Amax(rll ti)f 

S = n R 2 — 77J7T-. From (1571 ) and (|59l ) we note that we have the same minimum bit energy in the 
presence of QoS limitations while the wideband slope diminishes with increasing 9. 

When we have equal power allocation, i.e., = ^1, it can be immediately seen from the result of 
Corollary Q] that 

E b n T log e 2 



AWin E{tr(HtH)} 



(61) 



c = 2E 2 {tr(HtH)} 

° n if E{tr((HtH) 2 )} + fff(E{tr 2 (HtH)}-E2{tr(HtH)})' 

5 We note that the expressions in J54t and J56t differ from those in [2] by a constant factor due to our assumption that the units of Ce 
is bits/s/Hz/dimension rather than nats/s/Hz/dimension. 
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Assume that H has independent zero-mean unit-variance complex Gaussian random entries. Under this 
assumption, we have [3] 

E{tr (H+H)} = n R n T , E{tr 2 (HtH)} = n R n T (n R n T + 1), E{tr ((H t H) 2 )} = n R n T (n R + n T ). 

(63) 

Using these facts, we have the following minimum bit energy and wideband slope expressions for the 
uniform power allocation case when the entries of H are i.i.d. zero-mean unit- variance Gaussian random 
variables: 

Efj log 2 2 

— = — — and S = nR+nT eTB bits/s/Hz/(3 dB)/receive antenna. (64) 

No min n R R nT — + nT log ^ - 

We note that while the minimum bit energy depends only on the number of receive antennas, the 
wideband slope is a function of both the receive and transmit antennas. Note that the wideband slope 
expression is per receive antenna. Without this normalization, we have 

So = nR+nT 2 QT B— bits/s/Hz/(3 dB). (65) 

n r riT urut log e 2 

From (|65l) . we identify the interactions between the spatial dimensions and QoS constraints. Note that 
more strict QoS constraints and hence higher values of 9 tend to diminish the wideband slope. On the 
other hand, we see in the second term in the denominator of (1631) that the impact of the presence of 
QoS constraints is being diminished by the product of the number of transmit and receive antennas, 
n R riT- Hence, increasing the number of transmit and/or receive antennas can offset the performance 
loss due to queueing constraints. 

V. Minimum Bit Energy in the Wideband Regime 

In the previous section, we have assumed that the bandwidth of the system is fixed as the transmission 
power P diminishes and system operates in the low-power regime. Here, we study the regime in 
which the bandwidth increases while P is kept fixed. Note that as the bandwidth grows, the flat-fading 
assumption will no longer hold and the input-output relation given in (OQ) will not be an accurate 
description. On the other hand, if we decompose the wideband channel into parallel, noninteracting, 
narrowband subchannels each with bandwidth that is equal to the coherence bandwidth B c , then we can 
assume that each subchannel experiences independent flat fading and has an input-output relation that 
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can be expressed as 



yi = HjXi + ni i = l,2,...,m (66) 

where x, and are the input and output vectors of the i th subchannel, and Hj is the i th subchannel 
matrix, rij represents the additive zero mean Gaussian noise vector with ^{njiij} = N I in the i th 
subchannel. We assume that the input in the i th subchannel is subject to £'{||x i || 2 } < ^- where Pi is the 
power allocated to the i th subchannel. We assume that the number of subchannels is m and hence we 
have B = mB c and YlT=i P% = P where B and P denote the total bandwidth and power, respectively, 
of the wideband system. Under these assumptions, the maximum instantaneous transmission rate in the 
i th subchannel with covariance matrix K x i is 

B c log 2 det (l + TiRSNRiHiK^HjJ bits/s (67) 

where SNRj = ^ N ■ Due to the independence of fading in different subchannels, the total transmission 
rate over the wideband channel is 

m 

B c log 2 det + n K SNR i H i K a . )i HtJ bits/s (68) 

i=i 

which is achieved by independent signaling over different subchannels, i.e., by choosing {xj}™ 1 as zero- 
mean independent Gaussian vectors with covariance matrices {K x ^j™ v Then, for the transmission rate 
in (|68l) . the effective rate is given by 



1 i ( m 

Re(snr) = - eTBnR log e E jexp i-9TB c ^2^og 2 det (l + nijSNRjHjK^H 




(69) 



= -^^-log e fjE{exp (-0T£ c log 2 det (i + r^SNR^K,, t H{)) } (70) 
= -^Tg^^log e E{exp (-0TB c log 2 det (i + n^NR^K^Hl)) } (71) 

where (1701) follows from our assumption that {Hj} are independent subchannel matrices and the fact 
that the expected value of a product of independent random variables is equal to the product of the 
expected values of the individual random variables. In general, effective capacity can be obtained by 
maximizing the effective rate expression in (fTTT) over all power allocations {Pi} and covariance matrices 
{K^j}. If the channel is known at the transmitter, {Pi} and {K X| J can depend on the realizations of 
the channel matrices {H, }. 
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We simplify the above setting by assuming that HjK^jH] has the same distribution for all i = 
1,2, ... ,m. For instance, this assumption would hold when {Hj} are identically distributed, and K. Xii is 
the same fixed matrix for all i or is a random matrix with a common distribution for all i (e.g., K. x ti = 
uu', where u is the random eigenvector that corresponds to A max (H|Hj), has the same distribution for 
all i when {H;} are identically distributed). Under this assumption, we can eliminate the dependence 
of HjK^jHj on the time index i, and show from the concavity of the expression (I7TI ) with respect to 
signal-to-noise ratidj that the effective rate is maximized by having SNRj = n ^^ B = nR N B = SNR ^ or 
all i, i.e., by distributing the total power equally over the subchannels. Now, the effective rate expression 
becomes 

Msnr) = - 0T g rn log e E {exp (-9TB C log 2 det (i + t^snrHK^) ) } (72) 
= -— ± log e E{exp (-9TB c \og 2 det (l + nRSNRHK^H^) ) } (73) 

where we have used the relation B = mB c . 

Now, we analyze the effective capacity and energy efficiency in the wideband limit in three scenarios: 

1 ) Rich Multipath Fading: In a system with bandwidth B, the maximum number of resolvable paths 
is proportional to BT m = where T m denotes the delay spread and B c = In rich multipath fading, 
the assumption is that the number of independent resolvable paths increases linearly with increasing 
bandwidth. Therefore, in rich multipath fading, coherence bandwidth B c remains fixed as B increases 
while snr = -g^- diminishes to zero. Then, from the similarity of the effective rate expressions in (fT6l) 
and (1731) and the fact that B is fixed in (fT6l) in the low-power regime analysis, we immediately conclude 
that the wideband and low-power results are identical in rich multipath fading under the assumptions 
that lead to the effective rate expression in (1731) . 

2) Sparse Multipath Fading: In sparse multipath fading, it is assumed that the number of independent 
resolvable paths increases at most sublinearly with bandwidth [23] [24]. Hence, in this case, B c increases 
with increasing bandwidth. In the special case in which the number of resolvable paths is bounded, 
B c increases linearly with B while the number of subchannels m remains fixed. For instance, such a 
scenario is considered in [25]. For this case, we have the following result on the minimum bit energy 
required in the wideband regime. 

6 Since -0TB c log 2 det ^1 + nflSNRiHiKa^Ht J is a convex function of SNR for given HiK^Hj, 

-9TB C log 2 dotfl+nRSNRiHiK;^ ,Ht] . „ • -w , j 

e v is a log-convex function. Moreover, since log-convexity is preserved un- 

der sums [32, Section 3.5.2], E |exp ^—9TB C log 2 det ^1 + nRSNRiHiK^iH]^ | is log-convex, implying that 

log e E |exp (—6TB C log 2 det (l + nflSNRiHiK^iH]^ | is a convex function of SNR. Since the sum of convex functions is 
convex [32], and the negative of a convex function is concave, we conclude that the expression in d71b is a concave function of SNR. 
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Theorem 3: Assume that the number of independent resolvable paths remain bounded and fixed in 
the wideband regime as B increases. In this case, the minimum bit energy for a given co variance matrix 
Ka- is given by 

w — 

-^b _ mNp 

N Vmm _l oge ]E | e -S^k^2 tr ( HK - Ht )| 

When the channel is perfectly known at the transmitter, information can be sent in the maximal- 
eigenvalue eigenspace of H^H and the required minimum bit energy becomes 

OTP 

^b mN (15) 



No mm _ l oge ]E | e "S^toi^2 Amax ( HtH )| 



If only statistical information of the channel is available at the transmitter, the minimum bit energy can 
be obtained by minimizing (|74|) over all permissible covariance matrices, i.e., 

eTP 



5™ P — otp i . /ttt/ TTt ^ • (76) 



tr (K,)<i-log e E{e-^^ tr(HK ^ Ht) }' 
Proof: For a given input covariance matrix K. x , the bit energy required for reliable communications 



under QoS constraints is 

„ p 

&b _ SNR n R BN 



N R b (snr) -^^log e E{exp(-0T5 c log 2 det(I + n /? SNRHK :E Ht))} 

8TP 

mN 



- log e E jexp (-6TB C log 2 det (i + 



HK T Ht 



9TP 

mNo 



log e E{exp (-0T5 c £>g 2 (l + 



A (IIK II 



(77) 
(78) 

(79) 



mB c N " 

where Ai(HK a .H^) denotes the i th eigenvalue of the matrix HK^-H*. Above, (|78l) is obtained by using 
the relation B = mB c and performing some straightforward algebraic operations, and (1791 follows from 
the fact that det (A) = \\ i Aj(A). Note that under the assumption of fixed number of resolvable paths, 
B c increases linearly with B while m is fixed. Hence, only the denominator of (1791 varies with B. 
From the fact that the function xlog 2 (l + -) is a monotonically increasing function of x > for any 
constant a > 0, we can easily see that the minimum bit energy is achieved as B — > oo. Since B c also 
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grows without bound as B increases, we have 



E h 



v = lim r 7 ^7 xVT ( 8 °) 

N 0min B^oo _ ^ E | exp (_q TBc ^ ^ + _£_ A . (HKa;H t) j j J 



6>TP 
mNo 



log e E{exp (-^^E^A^HK^Ht))} 

6>TP 

m./Vo 

logeEjexp^g^Ei^CHK.Ht))} 



(81) 
(82) 
(83) 



- l°g e E {exp (-g^ tr (HK«Ht)) } ' 

(ED is obtained using the fact that as B c -> oo, we have £ c log 2 + Xi(HK x W) j -> j^^A^HKj: 

(f8~3l) follows from the property that J2i = tr (A). Note that (f83l) proves (|74l) which is the minimum 

bit energy for a given covariance matrix K x . 

Recall that it is shown in the proof of Theorem \T\ that 

tr (HK^H^) < A max (H t H) (84) 

and this upper bound can be achieved by transmitting in the maximal-eigenvalue eigenspace of Htf, 
e.g., by having K^. = iiu* where u is the eigenvector that corresponds to A max (H^H). If the transmitter 
perfectly knows the realizations of the channel matrix H, then this transmission strategy can be employed 
and the minimum bit energy becomes 

F ^ 

TT = 7 NT- (85) 

No - -logeElexp^g^A^HHt))} 

If the transmitter has only statistical knowledge of the channel matrix, the minimum bit energy can be 
determined by finding the input covariance matrix that minimizes (1831) . ■ 
Remark 3: By applying the Jensen's inequality, we can easily see that 

log p E( e -^^ tr{HK * Ht) ) > E(log e e"^^ tr(HKa;Ht) ) =e{- 9TP 1 trfHKgHt)! 
I J I J { mN log e 2 J 

(86) 

which implies that 

TP ^If- 1 O 

f± ^ mNo > lQ ge^ ( n 7) 

iVn min _l„„ W L-Si^»(HK I Ht)l -tr(HW 



^Omin _ log e E je - ^ 1 ^ tr ( HK - Ht )| 
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Similarly, we can show 



OTP 

> - "f T e " • (88) 



LP g£-£- 1 o 

-^6 _ mNo > [ °Se Z 

^Omin _ l g e ]E | e "S^T3i^2 Amax ( HtR )| ~~ A max (HtH) 



E b . ;^Vo . log e 2 log e 2 

iV 0min K?>rO 
ti(K x )<l 



^ mm = 

log e E {e-^^ tr(HK ^ Ht) } " ^gg, tr(HK«Ht) A max (E{HtH}) 



Note that the right-hand sides of the above inequalities are the minimum bit energy expressions in the 
low-power regime and also the wideband regime with rich multipath fading due to the equivalence of the 
two. From this, we immediately conclude that the sparse multipath fading with bounded number of re- 
solvable paths (or equivalently bounded number of subchannels) induces additional energy requirements 
in the presence of QoS constraints. 

Remark 4: Recall from the result of Theorem [2] that when the transmitter has only statistical knowl- 
edge of the channel, the optimal transmission strategy in the low-power regime (and also in the wideband 
regime with rich multipath fading) is to transmit the information in the maximal-eigenvalue eigenspace 
of E{H^H}. On the other hand, we note from Theorem [3] that this is not necessarily the optimal 
transmission technique in the wideband regime with sparse fading. The optimal input covariance is the 
one that minimizes (|74|) . Note further that for small we have the following first-order Taylor series 
expansion of the denominator of (|74|) : 



_ logeE{e -^.<H K ,„ t)} ^_l_ tt(HKiHt)+o Q. m 

Hence, when 9 or P is small or m is large, the input covariance that is optimal to the first order is the 
one that maximizes tr(HK x .H^), i.e., in this case, transmission in the maximal-eigenvalue eigenspace 
of E{H^H} is optimal as in the low-power regime. 

Theorem [3] holds for the case in which the number of resolvable multipath components remains 
bounded. Another scenario in sparse multipath fading is the one in which the number of resolvable 
paths increases with bandwidth but only sublinearly. In this case, both B c and m increase without 
bound as B — ► oo due to the sublinear growth of B c . Therefore, the minimum bit energy results can 
be obtained by letting m — » oo in the results of Theorem 

Theorem 4: Assume a sparse multipath fading scenario in which the number of independent resolvable 
paths increase sublinearly with bandwidth. In this case, the minimum bit energy for a given input 
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covariance matrix is given by 

E> = log. 2 
iV 0mln trfHK^Ht) 

When the transmitter perfectly knows the channel matrix H and when it knows only E{H^H}, the 
minimum bit energies are 

^ = lo ^ 2 and ^ = ^ (92) 

iVomin A max (HtH) N 0min A max (E{HtH})' 

respectively. 

Proof: As mentioned above, proof follows by finding the limiting values of the minimum bit energy 
expressions in Theorem [3] as m — > oo. For the case of fixed covariance matrix K. x , we have 

t± = ii m (93) 

-^Omin rn^oo TC I 17^72 te ( HK ^ Ht ) 



log e E | e ~ ^ ^ tr (HKt Ht ) } 



8TP 



lim ^ (94) 

^ W 1 tr(HK x Ht)+o(^) 



mTVo log e 2 



lim -, r (95) 



i ^tr(HK a; Ht) + ^ 



1 (96) 



8TP 

tr (HK,Ht) + lim. 



ttiNq 

log e 2 



tr(HK T Ht)- 



(97) 



(|94l ) is obtained by using the first-order Taylor expansion in (|90l ). (1931 ) follows by dividing the numerator 
and denominator by Finally, (|971) is obtained immediately from the definition that lim^o — 0. 
The expressions in (|92|) are determined as in the proofs of Theorems Q] and |2] by choosing the input 
covariance matrix as K. x = uu* where u is the eigenvector that corresponds to the maximum eigenvalue 
of H^H (when H is perfectly known at the transmitter) or of E{H^H} (when only E{H^H} is known 
at the transmitter). ■ 
Remark 5: Theorem |4] shows that as long as the number of subchannels m grows without bound in 
the wideband regime, the minimum bit energy requirements are the same as those in the low-power 
regime and wideband regime with rich multipath fading in which m increases linearly with bandwidth. 
Note that since each subchannel experiences independent fading, m can be seen as a measure of the 
degrees of freedom in the system. Therefore, if m is bounded, the degrees of freedom is also bounded 
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and that results in increased energy requirements as discussed in Remark [3j On the other hand, if the 
degrees of freedom increase with bandwidth, we have the same minimum bit energy values even though 
the increase is sublinear. However, for this case, we will observe in the numerical results in Section 
IVIII that approaching the minimum bit energy is very slow and demanding in bandwidth due to zero 
wideband slope. 

Remark 6: Note that having m — > oo for fixed 9 > in the minimum bit energy expressions in 
(|74|)-(|76|) is the same as letting 9 — ► for fixed m. Hence, even if m is bounded, the minimum bit 
energies given in Theorem @] are attained when = 0. This indicates that multipath spars ity does not 
affect the performance in the absence of QoS constraints. 

VI. The Impact of QoS Constraints in the High-SNR Regime 

In this section, we consider a single flat-fading channel and analyze how QoS limitations affect the 
performance in the high-SNR regime. In contrast to the previous sections where general models are 
used, we here consider a specific fading scenario in which the components of H are independent and 
identically distributed (i.i.d.) Gaussian random variables with zero mean and unit variance. Moreover, 
we assume that the input covariance matrix is K x = ^-1. Note that this covariance matrix is optimal 
in the sense of achieving the ergodic Shannon capacity when H has the above distribution and the 
transmitter does not know the realizations of H [26]. 

Now, for the considered channel and input models, the effective rate is given by 

#£,id(sNR) = -0-7g log e E jexp (^-9TB log 2 det ^1 + ^snrHH 1 "^ j J bits/s/Hz. (98) 

Note that in the above formulation, we have not normalized the effective rate expression with the number 
of receive antennas ur, and we have chosen a slightly different font from before and use the notation 
Rem to denote this unnormalized effective rate. 

As also pointed before, the effective capacity and effective rate expressions are proportional to the 
logarithm of the moment generating functions of instantaneous transmission rates. For the channel and 
input models considered in this section, Wang and Giannakis in [28, Theorem 1] provided an expression 
for the moment generating function of instantaneous mutual information. Applying this result to our 
setting, we obtain 

where T(-) is the Gamma function, k = min(nR, tit), and d = max(n#, nr) — min(nR,n;r). Moreover, 
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G is a k x k Hankel matrix whose (i, j) th component is 



oo / \ —0TB log 2 e 

nR 



9i,j = J [ 1 + ~ SNRZ J z i+j+d e- z dz i,j = 0,1,... ,k-l. (100) 

Therefore, we have 

„ , x 1 , / det(G(0,SNR)) . 

i?Eid(SNR) = --—log, ^ . (101) 

In order to quantify the impact of the QoS constraints on the performance in the high-SNR regime, we 
consider two measures, «Soo and C^, which are defined as 

i? g ,id(SNR) 

= hm — (102) 

SNR^oo log 2 SNR 



lin, llog 2 SNR- ^ (SNR) ). (103) 



and 



*^oo I "oz r* 

SNR^oo \ S, 

Note that while <Soo denotes the high-SNR slope in bits/s/Hz/(3dB), Cqo represents the power offset 
with respect to a reference channel having the same high-SNR slope but with unfaded and orthogonal 
dimensions [4]. With these quantities, the effective rate is approximated at high snrs as 

Rem = Soo(log 2 snr - C^) + o(l). (104) 

The above high-SNR approximation was first introduced and used in [27] in the study of code-division 
multiple access systems with random spreading, and was later employed in [4] in the study of ergodic 
Shannon capacity of multiple-antenna systems. Here, we apply this approximation to the multiple- 
antenna systems operating under statistical queueing constraints. The next result identifies the values of 
and for a subset of values of the QoS exponent 9. 
Theorem 5: Assume that the components of channel matrix H are independent and identically dis- 
tributed (i.i.d.) Gaussian random variables with zero mean and unit variance. If the QoS exponent 
satisfies 

max(B fi , n T ) - mm(n R , n T ) + 1 
TB log 2 e 

then, we have 

Soo = min(n R ,n T ), (106) 
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and 



lo S 2 £ + grfc lo Se ^ { e_9TB 1082 dCt HHt } "r < n T 



(107) 

riR > n T 

Proof: Note that the components of the Hankel matrix G, which appears in the effective rate 
expression in (11011) , can be written as 

roo / i „ \ -0TB log 2 e 

g iyj = s^R- eTBlo ^ e ( — + —^) z i+j+d e~ z dz i,j = 0,l,...,k-l (108) 



SNR 



tit 



where k = mm(riR, rvp). As snr — > oo, the integral in the above expression goes to a nonzero and finite 
value if —9TB log 2 e + i+ j + d> —1 since < J °° z a e~ z dz < oo for a > — 1 and J °° z a e~ z dz = oo 
for a < —1. Note that this condition is satisfied for alH, j — 0, 1, . . . , k — 1 by our assumption in (| 1051) . 
Now, we can immediately see that g it j for all i,j scales as SNR~ eTBl ° S2 e as snr — > oo. Therefore, the 
determinant of G scales as sNR~ fc6,TBlog2 e . This lets us conclude that 



R 



E,id 



SNR 



9TB 



1 , / det(G(0,SNR)) 
log e 



1 



n^r(d + <) 



log e (sNR- fceTBlos 2 e ) + 0(1) (109) 



9TB 

k(\og 2 e) log e snr + 0(1) 

k\og 2 snr + 0(1) 

min(nR, n T ) log 2 snr + 0(1), 



(110) 

(HI) 
(112) 



establishing that Soo = mm(nR,riT) for the values of 9 specified in the theorem. Above, 0(1) denotes 
the terms that approach a finite constant as snr — ■> oo. 

Next, we consider the power offset C^. Assume that n R < n T . Under this assumption, we have 

#£,id(SNR) 



SNR- 



lim logn SNR — 



SNR- 



lim 

SNR^oo 



lim 

SNR^oo 



lim ( loffn SNR — 



#E,id(SNR) 



( 



l 

9TB 



log 2 SNR + 



. _ f -0T J Blog,det('l+^B-SNRHHt') 1 \ 

log e E < e 2 v «t / ^ 



log 2 SNR + 



gj B log e E J" e - eT ^ 1 °g2SNR-erBiog 2 det( g ^ R i+^HHt) j > 



n R 



\ 



(113) 
(114) 

(115) 



(116) 
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/ 



lim 

SNR^oo 



-n R log 2 SNR + ^ log e E <J e 



log 2 SNR 



-erBlog 2 det(gjl E I+^HHt) \\ 



riR 



\ 



lim 1 log e E( e - erB1 ° S2dGt (sNR I+ ^ HHt ^ 
snr^oo 6TBn R Be \ 



1 



BTBn R 



, . —9TB logo det( — HH+ 

log e E<e 2 V"T / 



(117) 

(118) 
(119) 
(120) 



Above, while (11161) is obtained by noting that 

6TB\og 2 det (I + — SNRHH f ] = 9TB log 2 snr Hr + 9TB log 2 det ( — I + — HH f 
\ n T J \SNR n T 

the remaining steps follow through straightforward algebraic operations. The result for the case in which 
fix < n R can be readily proved by applying the above procedure to 

/ 



lim 

SNR^c 



log 2 SNR + 



1 l iu> f -eTBlog,det(l+^8.SNRHtH) 1 \ 

^log e Eje nr ^ J > 



\ 



(121) 



Remark 7: When ergodic Shannon rate (or equivalently effective rate with 9 = 0) is considered, it 
is well-known that the high-SNR slope for the i.i.d. Rayleigh fading channel is = mm(n R ,n T ). 
The above result shows that the high-SNR slope does not get affected by the queueing constraints when 



< 



maji(nfi,nj')- min(n^ ,nr)+l 
TBlog 9 e • 



Remark 8: For the case of 9 = 0, it is shown in [4, Appendix B] that the power offset in the i.i.d. 
Rayleigh fading is^ 



log 2 ^-^£{log 2 detHHt} n R <n T 
log.^-^-Ejlog.detHtH} n R >n T 



nR n R 
nj_ _1_ 

riR n T ' 

By Jensen's inequality and strict concavity of the logarithm function, we have 



(122) 



— log K/ e -^Blog 2 detHHt| > _J_ E (! -eT B log 2 detHHt| 

BriR I ) 9TBn R I J 

= — -£ {log 2 det HHH, for 9>0 

n R 



9TBn R 



(123) 
(124) 



7 In [4], signal-to-noise ratio is denned as SNR = ^frjr^ipn - ■ Due to the presence of ur in the numerator in the SNR definition, the 



first term of C^o in [4] is log 2 tit instead of log 2 — . 
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which shows from the comparison of (11071) and (11221) that the presence of queueing constraints result 
in higher power offset values in the high-SNR regime. 
Remark 9: Note that by Holder's inequality, we have 

{E{\x\ r }) l ' r < {E{\x\ s }f s (125) 



for < r < s. Note further that the second term in the expression of can be expressed as 







— \— log e E { e~ eTB det HHt ) = — log e ( E { e- 0TB dct HRt ) ) "* . (1 26) 



Application of the inequality in (11251) to \ E | e -^Bio g2 detHHt | j 0TB s jj OWS t ^ at power offset 
in a non-decreasing function of the QoS exponent 9. 

Theorem [5] characterizes and for a certain range of values of 9. The next result gives a partial 
answer to what is expected when 9 > max ( n R' n T)-^(nR,n T )+i ^ ^ considering the case of single-antenna 
transmission and reception, i.e., tit = % = 1. 

Theorem 6: In a Rayleigh fading channel with single transmit antenna and single receive antenna 
(i.e., n T = n R = 1), the high-SNR slope is 

Soo = — -J (127) 

fc'Ti? log 2 e 

when 6 1 > 



TBlog 2 e- 

Proof: When we have n T = n R = 1, the effective rate expression is 



R E ( sm ) = -J^log^le-^M^W*)} (128) 
^log„E{e'«.('+ s »wr"""°'} (129) 



9TB 



^log e E{(l + SNR|/ i | 2 )- eTBlog2e } (130) 



= ~WB loge l ( 1 +^r Sz£e ~ Zdz ( l31 ) 

where (11311) follows from our Rayleigh fading assumption which implies that z = \h\ 2 has an exponential 
distribution. Note that this effective rate expression can also be immediately seen to be a special case 

8 Without loss of generality, we consider the case in which ur < nr. 
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of the expressions in (1 1001) and (1 1 1 b ■ Now, we prove the result through the following steps: 

R E (sm) 



lim 

SNR^oo log 2 SNR 



lim 

SNR^oc 



lim 

SNR^oo 



dz 



9TB 



SNR r°° 
SNR JO 



log 2 SNR 

f™(l + SNRz)- eTBk,S2e e- z dz 



lim 

SNR^oo 



log 2 SNR 

log e SNR - ^rg log e SNR J °° (1 + SNRz)" eTjB1 ° g2 6 e' z dz 



lim 



hs log e SNR 9TB 



0TB 



SNR^oo log 2 SNR 
1 

+ lim 



+ 



log 2 SNR 

log e [SNR J °° (1 + SNRz)- 0TB1 ° g * " 6~ z dz 



log 2 SNR 

flfB lQ Se SNR Jo I SNR + Z J 



(132) 
(133) 

(134) 

(135) 

(136) 

(137) 

(138) 

(139) 
(140) 

Above, (11351) is obtained by multiplying the integral inside the logarithm in the numerator by as 
shown in (|134l) . and by using the fact that the logarithm of the division is equal to the difference 
of the logarithms. (11361) follows by separately writing the fractions. (11371) is obtained by evaluat- 



9TB lof 
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9TB lof 
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6TB lof 
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QTB lof 


, r 2 e 



lim 



+ lim 



log 2 SNR 

^ log e [snr-^ 1 ^ e+1 esHR r {-QTB log 2 e + 1, ^) 



log 2 SNR 



0TB 



loge 



e SNR 



i_ r( -0TB log a e+1, g±= 



-BTB log 2 e+1 

SNR 



log 2 SNR 



ing the limit of the first fraction, and by expressing (1 + snrz 



-0TB\ 



in the second fraction as 



SNR -tfTBlog 2 e(_l_ + ^ 



-8TBlog 2 e 
V SNR 



. (fT38l) follows from the fact that [33, Equation 3.382.4] 

6>TBlog 2 e / 1 x 

e~ z dz = e snr r I -9TB log 2 e, — 

\ SNR / 



(141) 



where T(a,x) is the upper incomplete Gamma function. (11391) is obtained by rearranging the terms in 
the numerator of the fraction in the second term. Finally, (|140l) follows by realizing that the limiting 
expression in (11391 ) is equal to zero. This is noted from the fact that as snr — > oo, we have 



e snr — > l 

r(-«rsiog 2e + i,i) _ 



l -9TB log 2 e+1 

SNR 



log 2 e - 1 



(142) 
(143) 
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indicating that the numerator in the limiting expression in (11391) is approaching a finite value as snr 
increases while the denominator grows without bound. The limit in (11431) is due to the fact that 

T(a,x) — 1 



as x^O (144) 



when a < 0, which is satisfied in our setting from our assumption that 9TB log 2 e > 1. ■ 
Remark 10: Theorem |6] shows for the single-antenna case that when 9 > max ( n fi' n r^-mm(n fi ,ra T )+i _ 
TB l og - , the high-SNR slope is iSoo = gTB \ og - < niin(n fl ,n T ) = 1, and diminishes with increasing 9. 
Note that by Theorem [51 Soo = min(n^, nr) = 1 when 9 < TB y og - in the case of single antennas at 
the receiver and transmitter. 

Remark 11: For the multiple- antenna case, we have the following additional discussion. An expres- 
sion for the components of the Hankel matrix G is given by [28] 

/oo / \ -9TB\og 2 e 

M+^snr^J z i+j+d e' z dz i,j = Q,l,...,k-l (145) 

71 

~ r {9TB log 2 e) sin (vr (d + i + j - 9TB log 2 e)) X 

V SNR ) 3 T(l + d + i + ]) / \ 

^jrprh—r—. „ Tn , s -xF 1 (l + d + i + j,2 + d + i + j-9TBlog 2 e,—?-) 

F(2 + d + i+ j — 9TB log 2 e) \ n R SNR J 

f^H T(9TB \og 2 e) / x 

- K r{ J . v— 1F1 9TB log 2 e,—d — i — j + 9TB log 2 e, — £- 

r (—a — i — J + 9TB log 2 e) \ n R SNR J 

(146) 

where jFi denotes the confluent hypergeometric function and has the following series expansion [33] 

1 F 1 {a,b,z) -^ — -1 + -- + -^-^-+ b(b + 1)(b + 2) 3J+... (147) 

Note that the expression in (11461) is valid when #T£> log 2 e 7^ ±(d + z + j) for all i,j because of the 
presence of the sinousoid in the denominator of the first term and the fact that T(x) = 00 or —00 when 
x is a negative integer. Under this restriction, we can see (by also noting that 1 F 1 (a,6, 0) = 1) that 
the first term inside the square brackets in (11461) scales as snr -1- ^ - - 7 while the second term scales as 
snr -6 * 71 - 8 log2 e as snr — > 00. Note that d = max(%, n T )— min(n^, ut) and i, j = 0, 1, . . . , mm(nn, n T ) — 

'The limit in J144t can be obtained from the following facts: A definition of the upper incomplete Gamma function is given by 
[33, Equation 8. 3514] T(a,x) = x a e~ x ^(l, 1 + a; x) = x a e~ x e~ xt (l + t) a ~ 1 dt. From this definition, we can easily see that 
lim^o ^ = / °°(1 + tr-'dt = ^ for a < 0. 
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1. Therefore, when 

9TB log 2 e > 1 + d + 2(min(raR, tit) — 1) = max(nR, + min(nR, tot) — 1, (148) 

the first terms with snr- 1 ^ - * - - 7 will dictate the rate at which gi/s approach zero for all Hence, 
we have 

-l-d-i-j 



(^Snr) ' 3 T(l + d + i+j) 
9i,j ~ r (9TB log 2 e) sin (vr (d + i + j - 9TB log 2 e)) r (2 + d + i + j - 9TB log 2 e) ' ' 49 ' 



as snr — > oo. Note that the matrix G, whose components g iy j are equal to the right-hand side of (11491) , 
is still a Hankel matrix as the components depend on the indexes only through (i+j). If the determinant 
is nonzero, it can be easily verified that the determinant of G scales as 

det(G) ~ snr-^-i (2j - 1} = S NR (mm{nR ' nT)) . (150) 



For instance, 

/ 

det(G) = det 



aSNR 1 6SNR 2 CSNR 3 
&SNR" 2 CSNR -3 (iSNR -4 



~ SNR- (1+3+5) = SNR" 9 (151) 



y |^ CSNR 3 (iSNR 4 eSNR 5 j J 

for large snr as long as the constant a, b, c, d, and e are such that det(G) is nonzero. Finally, we have 
under the aforementioned conditions that 



1 



/ det (G (9, snr)) \ (mm(n Rl n T )) 2 



-Rsid(SNR) ~ ~^r- log p r — — ~ ^ — — logo SNR, (152) 

,K ' 9TB \Uti^(d + i)J 9TB \og 2 e S2 



indicating that 



(mm(n Rl n T )) 2 
9TB log 2 e 



when 9TB \og 2 e > max(%,nT) + min(nR, n^) — 1. Note that under this condition on 9, Soo = 
UU 0TB\og2^ < m ^ n ( nR ^ nT )- Note also that the above conclusion reduces to the result of Theorem [6] 
when tir — tit = 1. 

VII. Numerical Results 

In this section, we numerically illustrate the analytical results obtained in the previous sections. In 
order to treat the Iow-snr and high-SNR regimes jointly, we consider the i.i.d. Rayleigh fading channel in 
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Fig. 1. Effective rate Re vs. SNR in the single-antenna case (i.e., when tir = tit — 1) for different values of 9 = 8TB\og 2 e - 



which the components of the channel matrix H are i.i.d. zero-mean, unit-variance, circularly symmetric 
Gaussian random variables. We further assume that the input covariance matrix is K x = ^1, and the 
effective rate is given by 



#£,id(sNR) = ~^7g log e E jexp (^-9TB log 2 det ^1 + ^snrHH 1 ^ ^ J bits/s/Hz. (154) 

Under these assumptions, we can easily compute the effective rate by using the formulation in (11011) 
and performing integral computations. We note that the computations of the effective rate in the 
correlated fading case can be done using the expressions of the moment generating function of the 
mutual information of correlated MIMO Gaussian fading channels provided in [29]. Summary of such 
non-asymptotic results, along with asymptotic spectrum theorems, on random matrices is presented in 
[30]. 

Figure \T\ plots the effective rate -R^id as a function of snr in the single-antenna case (n R = n T = 1) 
for different values of 9 = 9TB log 2 e. It is assumed that T = 1 ms = s and B = 100kHz 
= 10 5 Hz. Note that when 9 = or equivalently = 0, there are no statistical queueing constraints and 
the effective capacity is equal to the ergodic Shannon capacity. In Fig. [Q we observe that the effective 
rate in general diminishes with increasingly more strict queueing constraints (or equivalently higher 
9 values). As expected, under more strict buffer constraints, lower arrival rates are supported, and as 
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Fig. 2. Effective rate Re vs. bit energy ^ in the single-antenna case (i.e., when ur = tit = 1) for different values of 6 = 9TB log 2 e. 

a result, lower departure rates are seen. On the other hand, as predicted by the low-SNR results of 
Section [WJ all rate curves have the same slope at snr = 0. Note that this slope is the one achieved 
in the absence of QoS constraints (i.e., when 9 = 0). Therefore, the impact of queueing constraints on 
the performance lessens at low snr values. An intuitive explanation of this observation is that as power 
decreases, arrival rates that can be supported by the system diminishes as well, which in turn decreases 
the effect of buffer violation constraints. Note also that as discussed in Section|Vj results similar to those 
in the low-power regime are obtained in the wideband regime if the channel experiences rich multipath 
fading. Therefore, another interpretation of the above observation is that QoS constraints have less 
impact on the performance as the bandwidth increases in rich multipath environments. This is due to 
the fact that the number of noninteracting subchannels and hence the number of degrees of freedom 
increases with increasing bandwidth, and the system has increasingly higher diversity to combat with 
buffer constraints. 

Fig. Q] confirms the analytical high-SNR results as well. As predicted by Theorem [5l the high-SNR 
slope is the same as that achieved in the absence of QoS constraints as long as 9 = 9TB \og 2 e < 1. 
On the other hand, as proved in Theorem [6l high-SNR slope is strictly less than 1 when 9 > 1. The 
difference in the rates of increase at high snrs is clearly seen in Fig. [Q 

In Fig. [2l we plot the effective rate as a function of the bit energy in the single- antenna case. 
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E b /N (dB) 

Fig. 3. Effective rate Re vs. bit energy for 9 = 9TB log 2 e = 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 8 when ur = 2 and nr = 5. 



Confirming the discussion in Section IIV-BI we immediately note that the minimum bit energy for all 
values of 9 is —1.59 dB, which is the fundamental limit in the absence of QoS limitations. This is a 
consequence of the fact that the effective rate curves as a function of snr have the same slope at zero 
snr. However, since the second derivatives of the effective rate at snr = decreases with increasing 9, 
we observe in Fig. [2] that we have smaller wideband slopes, S , for larger values of 9. Similarly as in 
Fig. [Q we observe smaller high-SNR slopes, S^, when 9 > 1. 

In Fig. [3j effective rate vs. bit energy curves are plotted under the assumption that the number of 
receive antennas is tir = 2 and the number of transmit antennas is nx = 5. We still assume that T = 1 ms 
and B = 100kHz. In the figure, the curves from the top to the bottom are for 9 = 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 8 
in this order We again immediately note that the same minimum bit energy is attained for all values 
of 9 while the wideband slopes So are smaller for larger values of the QoS exponent. In this case, the 
minimum bit energy is . = 101og 10 = —7.61 dB At high snr levels, we observe that, as 

shown in Theorem [51 when 9 = 9TB log 2 e < max(n fi , n T ) — mm{n R , n T ) + 1 = 4, is the same as 
that achieved when 9 = (i.e., when 9 = 0). For 9 > 4, we note the gradual decrease in the high-SNR 

10 Note that when 9 — 0, effective capacity becomes equal to the ergodic Shannon capacity. For this case, rate is computed using the 
formulation provided in [25, Theorem 2]. 

"As opposed to f64b where = loge 2 , we have = loB § 2 in the figure since we plot the effective rate in bits/s/Hz without 

iv min n R Jv min n R 

normalization with the number of receive antennas. 
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Fig. 4. Effective rate Re vs. bit energy for nr = 2, 3, 4, 8, 15 when tir — 2 and (9 = 6TB log 2 e = 1. 



slope. 

When we compare Figs. [2] and [31 we see that the rate curves are much closer to each other in Fig. 
[3l indicating the resilience provided by spatial diversity against queueing constraints. This is further 
illustrated in Fig. HI where effective rate vs. bit energy curves are plotted for different number of 
transmit antennas when n R = 2 and 9 = 1. In this figure, we observe that the wideband slope S 
increases with increasing number of transmit antennas for a given QoS exponent 9. Moreover, we note 
that improvements are provided at all snr levels when the number of antennas is increased in the system, 
again pointing to the benefits of spatial diversity. 

Heretofore, the discussions on the Iow-snr regime apply to the cases in which the transmit power is 
small or the bandwidth is large but in a rich multipath fading setting. In Section |Vj we have remarked 
that sparse multipath fading has considerable impact on the performance in the wideband regime. In 
order to numerically illustrate these results, we provide Figs. \5\ and [6] In Fig. [51 effective rate 
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log e E 



jexp (- 



■9TB C log 2 det I I + — snrHH^ 
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(155) 



is plotted as a function of the bit energy. Above in (|155l) . B c denotes the coherence bandwidth, and 



SNR 



n R mB c No 



where m is the number of noninteracting subchannels, each experiencing i.i.d. zero- 
mean, unit-variance Gaussian fading. In this figure, we have ur = ut = 2, and = 10 4 , T = 1ms. 



33 



-8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 

E b /N (dB) 



Fig. 5. Effective rate vs. bit energy ^ for 6 — 0, 0.1, 0.5, 1, 2 in sparse wideband fading channels, np> = n T = 2. The number of 
subchannels is m = 5. The coherence bandwidth B c increases with increasing bandwidth. 

We consider the setting in which the number of subchannels is bounded while the coherence bandwidth 
increases with increasing bandwidth. We assume m = 5 and plot the curves by varying B c from 10kHz to 
10MHz. Therefore, bandwidth increases from 50kHz to 50MHz. As predicted by the result of Theorem 
[3l the minimum bit energy depends on 9 and increases with increasing 9. We note that for relatively 
large values of 9, considerably higher bit energies are needed when compared with the case of 9 = 0. 

In Fig. [5j we have assumed that the number of subchannels and hence the number of degrees of 
freedom is bounded, and B c increases linearly with increasing bandwidth. We have seen that having 
bounded number of degrees of freedom induces substantial energy penalty especially if the queueing 
constraints are stringent. Another scenario in sparse multipath fading is the one in which B c increases 
but only sublinearly with B. In such a case, the number of subchannels m increases with B as well. 
In Theorem HI we have shown for this scenario that the same minimum bit energy as in the case of 
9 = can be attained. This is depicted in Fig. [6l In this figure, the parameters are the same as in Fig. 
[51 except we now assume that m increases from 5 to 100 as B c increases from 10kHz to 10MHz. We 
note that in all cases, the minimum bit energy of —7.61 dB is approached. However, it is interesting 
to observe that the wideband slopes are zero when 9 > 0, indicating that approaching the minimum bit 
energy is very demanding in terms of bandwidth in the presence of queueing constraints. 
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VIII. Conclusion 

In this paper, we have investigated the performance of MIMO wireless systems operating under statis- 
tical queueing (or QoS) constraints, which are formulated as limitations on buffer violation probabilities 
in the large-queue-length regime. We have employed effective capacity as the performance metric that 
provides the throughput under such constraints. We have studied the effective capacity in the low-power, 
wideband, and high-SNR regimes. In the low-power regime, we have obtained expressions for the first 
and second derivatives of the effective capacity at zero snr under various assumptions on the channel 
knowledge at the transmitter side. We have shown that while the first derivative does not depend on 
the QoS constraints, the second derivative diminishes as these constraints become more stringent. As 
a byproduct of these results, we have demonstrated that the minimum bit energy requirements in the 
presence of QoS constraints in the low-power regime are the same as those required in the absence 
of such constraints. However, the wideband slope is shown to significantly get affected by queueing 
constraints. 

Results derived in the low-power regime are proven to apply to the wideband regime in rich multipath 
fading environments. On the other hand, we have noted that sparse multipath fading induces energy 
penalty if the number of noninteracting subchannels remains bounded in the wideband regime. In this 
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case, the minimum bit energy is shown to depend on the QoS exponent 9. If the number of subchannels 
increase with bandwidth but only sublinearly, we have seen that the minimum bit energy required in 
the absence of buffer constraints can be attained, but we have demonstrated in the numerical results 
that approaching this level is very slow. 

Finally, we have investigated the performance in the high-SNR regime by determining the high-SNR 
slope and power offset values. In particular, we have shown that if the QoS exponent is less than a 
certain thereshold, the high-SNR slope of min(n^, nr) can be maintained. However, in this case, we 
have remarked that there is still a price to be paid in terms of the power offset when queueing 
limitations are present. For the single-antenna case, we have proven that increasing 9 beyond a threshold 
starts affecting the high-SNR slope Soo. In such a case, Soo is shown to diminish with increasing 9. We 
have discussed extensions of this result to the multiple-antenna scenarios, and illustrated them through 
numerical results. 
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